Non-perturbative results for the spectrum 
of surface-disordered waveguides 

N. M. Makarov and A. V. Moroz 
Institute for Radiophysics and Electronics, 12 Acad. Proskura St., Kharkov, 310085, Ukraine 



Abstract 

We calculated the spectrum of normal scalar waves in a planar waveguide with 
absolutely soft randomly rough boundaries beyond the perturbation theories in the 
roughness heights and slopes, basing on the exact boundary scattering potential. The 
spectrum is proved to be a nearly real non-analytic function of the dispersion £ 2 of 
the roughness heights (with square-root singularity) as Q 2 — » 0. The opposite case of 
large boundary defects is summarized. 
In finite systems a long-distance (waveguide) signal propagation is caused by multiple 

reflections of the signal from opposite lateral boundaries. If the boundaries are irregular, 

each of the reflections is accompanied by a non-coherent scattering of the travelling wave. 

The multiple successive scattering events lead to substantial dephasing and attenuation of 

the primary signal. 

As far as we know, this effect was first consistently treated in worksBi. That simple 
and physically clear approach was based on the perturbation theory in the squared r.m.s. 
height ( of the boundary roughness and therefore required small enough (, viz. (k z () 2 <C 1, 
((/R c ) 2 -C 1, k( 2 /R c <C 1. Here k z is the transverse (normal to the waveguiding direction) 
component of the wavevector k and R c is the mean length of the boundary defects. In more 
recent paperJH] the theory of wave scattering from statistically rough surfaces was extended 
to arbitrary values of the Rayleigh parameter (k z () 2 . However, the other two inequalities 
of the above set were still necessary, which made it impossible to deal with, e.g., steep 



roughness slopes (((/R c ) 2 ^ 1) and/or the shadowing effect!'! (k( 2 /R c > 1). 

In this Letter we put forward the approach which is non-perturbative in the roughness 
heights and slopes. It is based instead on the exploitation of the exact boundary scattering 
operator. Due to this fact, we managed to extend the waveguide theory up to the quite 
general conditions of weak scattering (|3|), which are much less restrictive than the above ap- 
proximations. The most impressive advantage of our method is that it leads to new physical 
results even for the region of small heights ( {{k() 2 <C 1), where the waveguide propagation 
is believed to be well studied. The most surprising result is a non-analytic (square-root) 
dependence of the waveguide spectrum on the dispersion ( 2 (i.e. on the coefficient of cor- 
relation) of the roughness heights. This non-analyticity can not be in principle derived 
perturbatively. It means that the exact randomly rough boundary can not be reduced to 
the smooth random-impedance one even for the arbitrarily small irregularities. 

We consider a 2D (planar) strip confined to a region of the xz-plane defined by £(x) < 
z < d, where £(x) is a Gaussian distributed random function of the longitudinal coordinate 
x characterized by the properties 

(^(x)) =0, (t(x)Z{xf)) = C 2 W(|x-x'|). (1) 

The angular brackets denote an average over the ensemble of realizations of the profile 
function The binary coefficient of correlation W(|x|) has the unit amplitude (W(0) = 
1) and the typical width of R c . The spatial distribution of a scalar wave field inside the 
strip is governed by the Helmholtz equation and the temporal dependence by the factor 
exp(— iut), k = uj/c. Both lateral boundaries z = £(x) and z = d are supposed to be 
absolutely soft, i.e. the field vanishes upon them. We seek the averaged Green's function 
(Q(x, x'\ z, z')) to this Dirichlet boundary value problem. 

The exact integral equation for Q(x, x'; z, z') can be obtained through the use of Green's 
theorem: 

/oo 
dx s dz s G (\x -x 8 \;z, z s )E(x s , z s )Q{x sl x'\ z s , z'), (2) 
-oo 
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where S(x s ,z s ) is the effective scattering operator, 

d dC,(x s ) d 



l(x s , z s ) = 5[z s - £(x s )) 



dz s dx* dx. 



(3) 



and Qo(\x — x'\\ z, z') is the Green's function for the ideal (smooth) strip with £(x) = 0. 

To perform averaging of Eq. (|2|) we apply the elegant technique derived in worksii. 
Although the averaging procedure is rather straightforward, it has few fine points which will 
be discussed elsewhereH. As the result, we Bud the averaged Greeu's fuuctiou: 

/OO (fa 
—^exp[ik x (x-x')}G(k x ;z,z r ), (4) 
-oo AH 

where longitudinal Fourier transform G(k x ; z, z') is given by 

K x] z " z) - i-k 2 cot(k 2 d)M(k x y 1 j 



Here k z = k z {k x ) = yk 2 — k%, and Go(k x ; z, z') is the Fourier transform (similar to Eq. (f|)) of 
Qq{\x — x'\; z, z'). The self-energy M(k x ) is defined via the binary correlator of the scattering 
operator (|]): 

sm(k z z s 



M(k x ) — I dz s dz' s dx s z s exp(—ik x x s 

J — OO rh-. 



x (5(3?,, z s )G Q (\x s - x' s \; z a , z' s )E(x' s , z' s )} exp(ik x x' s ) Sm ^ zZ * . (6) 

By equating the denominator of Eq. (|) to zero, we obtain the dispersion equation for 
the rough-bounded strip. The solution to this equation in the lowest (linear) order in M(k x ) 



is k x = k n + 5k n . Here k n = y k 2 — {itn/d) 2 is the unperturbed longitudinal wavenumber of 
an n-th propagating normal waveguide mode and Sk n being the complex modification to k n 
caused by the wave scattering from the lower irregular boundary, 

5k n = j n + i{2L n )~ l = -M(k n )(Tm/d) 2 /k n d. (7) 

The real part 7„ of 5k n is responsible for variation of the phase velocity (for dephasing), 
while L n has the meaning of the attenuation length for the n-th mode. 



3 



In fact, a form of the solution (|5]), (^) is common and well-knownBB. Our improvement 
lies in the self-energy M(k x ) (|6]). This expression is nothing else but a first non- vanishing 
(quadratic) term in an expansion of M(k x ) in powers of the exact scattering operator 




We stress that such an approximation is essentially different from and is much more general 
than an extensively exploited first term in an expansion of M(k x ) in powers of the dispersion 



To find the domain of validity for Eqs. (|6]), (|7|) one can use ideas proposed in the bookO. 
We have proved that this domain coincides with the two natural requirements of weak 
wave scattering from a rough boundary, which are equations themselves with respect to the 
external dimensionless parameters {k() 2 , kR c , kd/ir, and n, 



The first of these implies smallness of the complex phase shift over the distance A n = 
2k n d/(irn/d) passed by an n-th mode between two successive reflections from the rough 
boundary Under this condition the mode experiences a large number L n /A n ^> 1 of the 
reflections before its amplitude will be substantially attenuated. Also this ensures smallness 
of 5k n in comparison with the unperturbed wavenumber k n , because the inequality k n A n > 1 
always holds. 

The second of Eqs. @ indicates that the phase shift must remain small over the typical 
variation scale R c (effective correlation radius) of the boundary scattering potential. Obvi- 
ously, this is the necessary and sufficient condition for correctness of statistical averaging over 
The quantity R c is defined as the typical width of the correlator (E(x s , z s )E(x' s , z' s )) 
in Eq. as a function of x s — x' s . It does not generally coincide with the mean length R c 
of boundary defects. 

We next calculate the correlator and the integrals over z s and z' s in Eq. (^|). Then we 
substitute the result into Eq. (0) and extract real and imaginary parts of 5k n . Finally, we 
get explicit formulae for 7 n and L n : 



c 2 . 



\Sk n \A n < 1 



\5k n \R c < 1. 



(8) 
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In 



C (irn/d) 2 ^ {im'/d) 



E 



+2C 



2 k n d n / = \ k n 'd 
2 (7in/d) 2 ^ [im'/df 



W s (k n , k n ,) - W s (k n , -k n , 



{im/df 
kn.d 



M 2 {k n ) 



E 



knd n , =nd+1 \k n '\d JO 



dx exp(—\k n i\x)Re exp(— ik n x)W(k n , i\k n '\; x) 



(9) 



L -i = C 2 cw^! g (™'M 2 



^n<^ n '=l kn'd 



W c (k n , k n ') + W c (k n , -k n >) 



(10) 



Here integer = [kd/ ir] is the number of the propagating normal modes in the smooth strip. 
The function W(k x , q x \ x) is the generalized coefficient of correlation (W(k x , q x ; x) ~ W(|x|) 
as C 2 ^0), 



W(k x ,q x ;x) = (4k z q z ( 2 ) 



(k z + q z f + (k s + q z )(k x - q x ) ( - — j - (k x - q^ 2 ' ' W '' 



(** - q z f + {k z - q z )(k x - q x ) [ + — ) + (k x - q x 



k z q z 

> k x q x 



-\-2{k x q x 



k. 



kr q 
k x q x 



q x k x -\- ( k x q x t 
q z & z rC z q. 



k z q 

S{k z + q z , k z - q z ;x) 



S(k z + q z ,k z + q z ;x) 
S{k z — q z , k z — q z ; x) 



(11) 



S(ti,t 2 ;x) = (tit 2 ) _1 sinh [m 2 C 2 W(|x|)1 exp \-{tj + 1 2 2 )( 2 /2 



(12) 



where q z = k z (q x ) and the functions Ws(c){k x ,q x ) stand for sine and cosine Fourier trans- 
forms of W(k x , q x ; x) respectively. The component M 2 (k x ) of the self-energy is given by 

k 2 00 



M 2 (k 



^2S(k z + nri' /d, k z — nn' /d; 0) 



2k 2 d n , =1 

S(k z + 7m' /d, k z + nn' /d; 0) — S(k z — nn'/d, k z — 7m' /d; 0) 
k. Tin 1 



k 2 d 



[S(k z + 7m' /d, k z + nri/d) 0) - S(k z - 7m' /d, k z - 7m' /d; 0)]}. 



(13) 



An essential distinction between 7„ @ and L^ 1 ([10|) is that the latter is formed by 
scattering of a given ra-th propagating mode into propagating waveguide modes with n' < 
only, while the former has much more complicated structure due to contributions of both 
propagating and evanescent (n' > rid) modes. This feature is a basis for surprising properties 
of 7n- 



Brief analysis and discussions: We start with a relatively simple and widely used limiting 
case of small boundary perturbations, when 

(k() 2 « 1. (14) 

Here Eq. (Jjl]) for L n is simply reduced to the standard result from the earlier wor k M by 
replacing W(k x ,q x ;x) with W(|x|). On the contrary, j n (§) shows an unconventional type 
of behavior. The reason is that the last two terms of Eq. (|9|) are mainly formed by those 
evanescent modes whose normal wavelengths (im'/d)' 1 are of the order of the roughness 
height (. Each of such 'resonant' modes gives a contribution to 7„ proportional to £ 2 , while 
the number n' of those modes is ~ d/( 3> n^, i.e. it is inversely proportional to (. All this 
gives a linear dependence of 7„ @ on the roughness height (, 

7 n ~ C {^n/df/k n d. (15) 

This formula is the main result of the Letter. It leads to the following significant conclusions: 

1. Since L^ 1 oc ( 2 as ( 2 — > 0, then 7 n L^ 1 and, hence, the entire spectrum shift Sk n 
(0) turns out to be nearly real, i.e. 5k n ~ 7 n oc (. This means that a signal propagating 
through a nearly smooth waveguide is dephased (chaotized) much earlier (over much shorter 
distances) than its initial amplitude is damped. 

2. From Eq. (|i~5|) and item 1 it follows that Sk n is a non-analytic (square-root) function 
of the dispersion £ 2 , or of the binary correlator ([]]): 5k n oc (C 2 ) 1 ^ 2 - 

3. Usually it is believed that the condition (|14D is sufficient to infer that any long- 
wave normal mode with {k z Q 2 1 (i.e. with (jm/d)^ 1 ^> () is mainly scattered into the 
long-wave modes as well, (irn'/d)^ 1 ^> This assumption allows immediate replacement 
of the exact Dirichlet boundary condition, formulated on a randomly rough boundary, by 
an approximate impedance one laid down on the averaged (deterministic) boundary z = 
with the random impedance ^(x). We have proved that such a reduction is groundless, 
because the 'resonant' evanescent modes with (irn'/d)^ 1 ~ ( dominate. Thus, the problem 
of wave propagation through a waveguide with an absolutely soft random boundary can not 
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be reduced to that with the smooth random-impedance boundary even for the arbitrarily 
weak perturbations. 

An important step in analyzing the case ( [H]) is to obtain the explicit weak-scattering 
condition. To this end, we substitute Eq. (|1|) into Eqs. (||) and apply the asymptotic 
R c ~ R c , which is valid for (k() 2 < 1. As a result, we find that Eqs. (§) can be rewritten as 

(fc^) 2 = (rcn/dfC 2 « min {l, (A n /R c ) 2 } . (16) 

This inequality is automatically satisfied within Eq. (|14[) if successive reflections of an n-th 
mode from the rough boundary are not correlateM^ (R c <C A n ). However, if the correlations 
are strong (A n <C R c ), then Eq. (|T6|) supplements Eq. (|14]) and may even become more 
restrictive than Eq. (|14]). Note that the roughness slope (/R c may far exceed unity within 
the limit (0), (|TJ). 

The situation with large boundary defects, when 

{Kf » 1, (17) 

is much more diverse and complicated than the case ( |H| ) . For the most part it is analyzable 
only numerically. Therefore we postpone the details until the longer paperS and list here 
only few of intriguing results: 

1. In contrast to the case (|14D , the imaginary part of 5k n (ff|) may well compete with its 
real part. Moreover, the situation with L^ 1 > \j n \ is rather typical. 

2. The real spectrum shift 7„ (^) reverses the sign as k( reaches some threshold value 
k( w 1.5 -T- 2.5, which is slightly dependent on the other parameters. So, the large boundary 
defects (|17|) may not only decrease, but also increase the phase velocity of a propagating 
wave. 

3. As k( < 2, 5k n is almost insensitive to the slope (/R c and to the presence of the 
shadowing effect, which is controlled by the Fresnel parameter k( 2 /R c . However, as k( 2, 
the weak-scattering conditions (§) are fulfilled for not too steep slopes only ((/R c ~ 2 3), 
while the strong shadowing effect is still allowed (k( 2 /R c ^ 8 4- 10). 



We mention that the approach presented in this Letter can be extended to vector wave 
fields. 

The authors thank Prof. A. A. Maradudin and Dr I. V. Yurkevich for useful comments 
on this work. 
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